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Abstract. Let p > 2, n\ < ■ ■ ■ < n p be positive integers and B\, . . . , B\ x ; . . . ; B\, . . . , B% be 
independent planar Brownian motions started uniformly on the boundary of the unit circle. 
We define a p-fold intersection exponent 5(711, . . . , Tip), as the exponential rate of decay of the 
probability that the packets (Jj=i Bj[0,t 2 ], i = l,...,p, have no joint intersection. The case 
p = 2 is well-known and, following two decades of numerical and mathematical activity, Lawler, 
Schramm and Werner (2001) rigorously identified precise values for these exponents. The 
exponents have not been investigated so far for p > 2. We present an extensive mathematical 
and numerical study, leading to an exact formula in the case n\ = 1, n<i = 2, and several 
interesting conjectures for other cases. 



1. Introduction 

1.1 Motivation and overview 

Finding exponents, which describe the decay of some probabilities, and dimensions of some sets as- 
sociated with stochastic models of physical systems is one of the core activities in statistical physics. 
While in general one often has to resort to numerical methods to get a handle on the values of the 
exponents, for planar models conformal invariance may help to answer these questions explicitly, and 
there is now a substantial body of rigorous and non-rigorous methods available. For example, by 
making the assumption that critical planar percolation behaves in a conformally invariant way in the 
scaling limit and using ideas involving conformal field theory, Cardy [Ca92] determined the asymptotic 
probability, as N — > 00, that there exists a two-dimensional critical percolation cluster crossing a rec- 
tangle. A rigorous proof of Cardy's formula was later given by Smirnov [SmOlj . Following considerable 
numerical work, see for example [LR781 IVo84] and references therein, Saleur and Duplantier [SD87] 
predicted the fractal dimension of the hull of a large percolation cluster using a non-rigorous Coulomb 
gas technique. Rigorous versions of this result have been given based on Cardy's formula, for example 
by Camia and Newman [CN061 ICN07] . 

In [DK88] Duplantier and Kwon suggested that ideas of conformal field theory can also be used to 
predict the probability of pairwise non-intersection between planar Brownian paths. Early research by 
Burdzy, Lawler and Polaski [BLP89] and Li and Sokal [LS90J was of numerical nature, but ten years 
later, Duplantier [Du98j gave a derivation based on non-rigorous methods of quantum gravity, and 
soon after that Lawler, Schramm and Werner [LSOlal ILSOlbl ILS02] gave a rigorous proof based on the 
Schramm-Loewner evolution (SLE), one of the greatest achievements in probability in recent years. 
We also mention here some very recent developments with the long term aim of making the quantum 
gravity approach rigorous, see Duplantier and Sheffield [DS08J, and Rhodes and Vargas [RV08J. 
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In this paper we look at joint intersections of three or more planar Brownian paths, a question 
which has been neglected so far in the literature, but which came up in our recent investigation of 
the multifractality of intersection local times |KM05] . In the simplest case, given three independent 
Brownian paths B , B 2 , B 3 started uniformly on the unit circle, we are interested in the asymptotic 
behaviour, as t — > oo, of the non-intersection probability 

F{B 1 [0,t]nB 2 [0,t]nB 3 [0,t] = 0}. 

Observe that this probability goes to zero, for t f oo, as three, or any finite number, of Brownian 
paths in the plane eventually intersect, see e.g. [MP 091 Chapter 9.1]. Recall for comparison, that the 
non-intersection exponents for three Brownian paths studied in the aforementioned papers deal with 
pairwise non-intersections, i.e. in the case of three Brownian motions either with 

F{B 1 [0 1 t]nB 2 [0,t] = 0, B 2 [0,t]nB 3 [0,t] =0, B l [0,t] nB 3 [0,t] =0}, or with 

P{B x [0,t] n (B 2 [0,t]U B 3 [0,t}) = 0}. 

Our study starts with the observation that, for positive integers rai,... ,n p and independent planar 
Brownian motions 

B\, ■ ■ ■ , B^; . . . ; Bf, . . . , B^ p , 

nontrivial exponents 

m n p 
? (n 1; ...,n p ) = - Urn — - logPj |J B}% t] n . . . n |J fl?[0, t] = 0) 

exist, see Theorem Q] and the subsequent remark. In Theorem [2] we show that, for 2 < 113 < ■ ■ ■ < n p , 
we have 

?(l,2,n 3 ,...,n p ) =2. 

These are the only exponents we could determine exactly beyond the well-known case of p = 2. 
Rigorous proofs of both theorems are given in Section [2j 

The bulk of this paper is devoted to the presentation of a detailed numerical study of the values 
of the, in our opinion, most interesting remaining exponents, see Section [3l One of the motiva- 
tions of this study was to test the conjecture, motivated by Theorem [21 that the value of the ex- 
ponents ?(nx, rt2, TI3, ■ ■ ■ , Tip) depend only on the two smallest parameters. This conjecture was not 
supported by our numerical investigations. 

Finally, we remark that we have not been able to use either SLE techniques or quantum gravity to 
derive even a non-rigorous exact prediction of the exponents if p > 2. We hope however that our 
numerical study triggers interest in this problem and that, as in the motivational examples discussed 
above, future research will address the question of exact formulas for multiple intersection exponents. 

1.2 Statement of the main theorems 

Let p > 2 and n\,...,n p be positive integers and B\, . . . , B\ ; . . . ; Bf, . . . , Bn p independent planar 
Brownian motions started uniformly on the unit circle dB(0, 1). We define p packets by 

<BV) := (jB}[0,r}(r)), ... ,&{r) := (J i?[0,Tf(r)] , 
where rj(r) := inf{t > 0: |Bj(t)| = r} and r > 1. 
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Theorem 1. The limit 

dni,...,n p ) := - lim — ^ — logPlosVr) n . . . n <B p (r) =0) 

r— »oc — log r y- J 

exists and is positive and finite. 
Remarks: 

• Using a standard argument, see [La96, Lemma 3.14], one can replace the paths stopped 
upon hitting the circle of radius r, by paths running for t = r 2 time units. This leads to the 
characterisation of the exponents given in the overview. 

• For p = 2 all exponents are known, see [LSOla, LSOlb, LS02]. The technique used to identify 
the exponents, which is based on the Schramm-Loewner evolution (SLE), does not seem to 
allow us to identify the exponents for p > 2. 

• We conjecture that one can strengthen this result, as this was done for p = 2 in [La95], and 
show that there exists a constant c > 0, depending on the starting points, such that 

lim r ^i>->™i>) p/BVr) n . . . D <B p (r) = j = c. 

However, this is quite subtle and would go beyond the scope of this paper. 

There is a trivial symmetry of the exponents, namely for every permutation a G Sym(p), we have 

?(ni, . . . ,n p ) = s(n CT (i), . . . ,n CT (p)). 
Moreover, there are two trivial monotonicity rules for these exponents 

(A) q(ni, ...,n p )< q(m, . . . ,n p _i), 

(B) . . . ,n p ) < ?(mi, . . . ,m p ), if raj < m* for i = 1, . . . ,p. 

As a result of the symmetry of the exponents, we may henceforth assume that the arguments of the 
exponents are increasing in size, i.e. m < • • • < n p . There is one interesting situation in which we can 
determine the exponents explicitly. 

Theorem 2. We have 2, n%, . . . , n p ) = 2 for any p > 2 and 2 < < ■ ■ ■ < n p . 
Note that to show this, by the monotonicity rules, it suffices to show that 
(1.1) ?(1,2, P ~- ,2) = 2. 

The proof of this fact is based on the technique of hitting the intersection of p — 1 Brownian paths by 
a further path, using an idea of Lawler, see [La89j or |La91j Section 3], originally used to determine 
the exponent <j(l, 2) = 2. 

1.3 Conjectures 

In this section we formulate the main conjecture motivated by our numerical studies. A detailed 
description of these studies and their outcomes will be given in Section [3l 

Let p €N and n±, . . . , n p EN with ri\ < ni < . . . < n p . Define 

k := min {£ £ {2, . . . , p} : n e+ i > n e } , 
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with k := p if the set is empty. We conjecture that 

(1.2) s(m,...,n p ) =q(n 1 ,...,n k ). 

In fact, this holds, by Theorem [2] for the case k = 2, = 2, and we have numerical evidence for 

• ?(1, 1, 2) = 1.2503 ± 0.0011 to be compared with <j(l, 1) = | 

• 1, 1, 2) = 1.02 ± 0.004 to be compared with ?(1, 1, 1) = 1.027 ± 0.005 

• ? (1, 3, 3) = 2.688 ± 0.01 to be compared with <r(l, 3) = 13 +v^ = 2.693000 . . . 

• ?(2, 2, 3) = 2.937 ± 0.01 to be compared with ?(2, 2) = f§ = 2.91666 . . . 

• ?(2, 3, 3) = 3.767 ± 0.06 to be compared with ?(2, 3) = 47+5 /^ = 3.738334113 . . . 



2. Proofs of Theorems Q] and [2j 

2.1 Proof of Theorem Q] 

Denote by x = (x\, . . . , x\ x \ . . . ; x\, . . . , xiJj vectors with ni -| + n p entries in M 2 , playing the role 

of configurations of our motions at time zero. Consider 

a r := sup P^V) H • • • fl <B p (r) = 0}, 

1*3 l=i 

where the subindex of P indicates the starting points of the Brownian motions. Using the strong 
Markov property and Brownian scaling, we get, for any r, s > 1, 

n\ Tip 



a rs < sup PJ \jB}[0,rj(r)]n---n\jBP[0,Tf(r)] 
[*$[=! L ,=l j=1 



ni ™p 

[j B}[r}{r),r}{rs)\ n ■ ■ ■ n (J ^[r?(r) )7 f («)] = 0} 



3=1 



3=1 



rti «p 

sup E* [l{ (J 5] [0, r) (r)] n • • ■ n (J [0, rj(r)] = 0} 
41=1 S=i i=i 

x P (BKrj(r)) { U B)[r}(r\r){rs)\ n • • • fl |J ^[tJ (r),7f(r.)] = 0} 



< a r a, 



Hence the function given by 6j := log a 2 * is subadditive and, by the subadditivity lemma, see e.g. [La91j 
Lemma 5.2.1], we thus have limj^oo bt/t = inf/,>o bt/t. Therefore, 

<f(ni,...,n„) := - lim —— log sup PjQSVr) n • • • n <B p (r) = 0} 
r-oologr lxi]=1 



exists, and is positive. 
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Next, we show that we can replace the optimised starting points by starting points uniformly chosen 
from the unit circle. Clearly, we have 

(2.1) P{«8 x (r) n • • • n <B p (r) = 0} < sup P^V) n • • • n <B p (r) = 0}, 

[*}[=! 

where P refers to the original scenario of Brownian motions started uniformly on the unit circle. 
Conversely, using the Markov property, for r > 2, we have 
sup PzjQSV) n • • • n <B p (r) = 0} 

< sup ix^^^f U B l[^ 2 )^i Wl " ' ' ' n U 5 i t 7 ! Wl = 

1-51=1 S=i J= i 

By the Harnack principle, the law of the vector (£?*-(-rj(2))) is bounded, uniformly in x, by a constant 
multiple of the uniform distribution on the (m + • • • + n p )-fold cartesian power of the circle dB(0, 2). 
Denoting this constant by C and using Brownian scaling, 

p { U B ) [°. T ) ( r / 2 )] n - n (J flj [0, (r/2)] = 0} 

(2.2) i =1 i =1 

>C _1 sup P x {«B 1 (r)n---n*B p (r) = 0}. 

Hi=i 

Combining and ([23]) yields that 

?(ni,...,np) := - lim —— logP{53 1 (r)n---n«B p (r) = 0) 

r^oo log V 

exists and coincides with q(nx, . . . ,n p ). Note, finally, that the monotonicity rule (A) implies that 
<r(ni, . . . , n p ) < s(ni, nz) < oo, and hence the exponents are positive and finite. o 



2.2 Proof of Theorem [2j 

Recall that it suffices to show (jl.ip . We start by formulating the key lemma. We let W , . . . , W p be 
independent Brownian paths. For r, s > denote by r l (x,r) the first hitting time by the motion W l 
of the circle dB(x, r) with centre x and radius r, and let t 1 (x, r, s) be the first hitting time of dB(x, s) 
after t 1 (x, r). 

Lemma 3. Fix x G £>(0, 1). Suppose that W 1 , . . . , VF P are independent Brownian paths started uni- 
formly on the circle dB(0,2). Define the set 

p 

(2.3) 2B:= P) W^'[0,r J '(0,4)] 

i=2 

and i/ie events 

E x , r = {W 1 [0,T 1 (x,r/2)]nW = (b}, 

(2.4) Ar x , r = {W 1 [O,T 1 (a;,r/2,r)]narr^0} ) 

ffx,r={ r W/2) <r i (0,4) /or a// i=l,...,p}. 
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Then 



lira inf ■ 

r|0 



logr| 



log¥[E x>r nN x>r \H x>r ] > 2, P ---, 2). 



Let us first see how follows from this lemma. Let 

r = inf{i >0 : W l {t) G2U}. 

Now let 25 be a collection of pairwise disjoint discs of fixed radius < r < 1/2 with centres in the disc 
13(0, 1), which has cardinality at least (2r)~ 2 . Then, obviously, 

i > p{w 1 [o,t 1 (o,4)] nau^0} > ^p{f 1 (T)eB,T<T 1 (o,4)}. 

Now, fix a disc B = B(x,r) G 25. The event {H /1 (r) £B,t< r 1 (0,4)} is implied by the events 

£^niVn {^0*^/2) <r 1 (0,4)}. 

Recall that 

F[H x , r ] =P{r 1 (x,r/2) <r 1 (0,4)} p = r°«. 
Combining this with Lemma El for any e > and sufficiently small r > 0, 

P{WV) G B, r < ^(0,4)} > r f(i.2,-.2)+s > 

This implies 

1 > ^ r ?(l,2,...,2)+ £ > r -2+?(l,2,...,2)+2 e ^ 

and therefore 2, . . . , 2) > 2 — 2e. The lower bound follows as e > was arbitrary, and the upper 
bound in (jl.ip follows from 2, . . . , 2) < 2) = 2, as is known from [La891 lLa9"T] . o 



Proof of Lemma El Before we describe the technical details we sketch the idea of the proof. Since 
the paths of p planar Brownian motions intersect with positive probability, by Brownian scaling, the 
conditional probability of N x r given H x ^ r is bounded from below as r — > 0. Hence this condition can be 
neglected when computing the probability in Lemma El For j = 1, ... ,p we decompose the paths W 3 
into the pieces W 3 [0,T 3 {x,r/2)] and W j [r j (x, r/2), t?(0, 4)]. By time reversal for W^[0,r^(x,r/2)], 
we can compare the probability in question with the non-intersection probability for packets of size 
ni = 1, ri2 = ■ ■ ■ = np = 2, which is of order ~ ^C 1 ' 2 '— > 2 ). 

We now come to the technical details, see the appendix in [MS09] for the necessary facts about Brow- 
nian excursions between concentric spheres. Let q 1 = r and gp = r/2 for j = 2, . . . ,p. Conditioned 
on {t 1 (x,q 3 /2) < t*(x,3)} the path W l [0, t 1 (x, q 3 /2)] is contained in an excursion from dB(x,3) to 
dB(x, $ /2). The time-reversal of this excursion is contained in the path of a Brownian motion W l 
started uniformly on dB(x, Q 3 /2) and stopped upon reaching dB(x,3), say at time r*(x,3). Analo- 
gously to and (f2T4|) define the set 

p 

W= p| (W j [0,T j (x,3)}UW j [T j (x,r/4,r/2),T j (0,4:)]), 

3=2 
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and the events 

E Xtr = {W 1 [0,T 1 (x,3)]nW = fH}, 

N x , r = { H W 3 [t 3 \xJ/2),t 3 {xJ/2J)] ± 0}, 

H x>r = {t 3 (x, q>/2) < r j (x, 3) for all j = 1, . . . ,p}. 

Note that W 1 [0,t 1 (x, p 1 )] nB(i,r/2) = and W j [Ti(x,pt /2),r i (x,p' /2,/jf)) C B(x,r/2) for j = 
2, . . . ,p. Hence 

p 

W 1 [0,T 1 (x,p 1 )]n(W\W) cW l [0,T 1 (x,p 1 )}nf]W j [r j (x,p i /2),T j (x,p j /2,p j )) = (D 

i=2 

which implies E x r C is^ r . Note that trivially, we have H xr C H x>r and iV^ r C N x>r which implies 

(2.5) E x>r Pi N Xji . Pi D -Ex.r H iV^r D -f/z,,-- 

Finally, note that 

_ Pffl'xr] _ Fjr^x,^ 1 ^) < t 1 (x,3)} p 1 
(2 " 6) /(X ' r) P[^J ~ P{ri(x,r/2)<ri(0,4)} p " 2 

for all x and for sufficiently small values of r > 0. 

By (|2.5p . (|2.6p and the definition of the conditional probability, we conclude 
(2.7) ¥[E Xtr n N x>r j H Xjr ] > f{x,r)F[E x>r n N x , r \H XyT \ . 

Fix e > 0. Invoking the definition of the exponent, the Harnack principle and Brownian scaling, for 
sufficiently small r > 0, 

V[E x>r \H x , r ] > P «iA-.a)+«. 

Define the compact sets 

C := {y = (y 1 , . . . ,y*>) : G dB(0, s> /2) for j = 1, . . . ,p} and 
D := {z = (z 1 ,...,z p ) : z 3 € dB(Q,Q>) for j = l,...,p}. 

For y G C and 2: G D let (W^, j = 1, . . . ,p) be an independent family of Brownian motions where 
each motion W 3 is started at y 3 and is conditioned to leave B(0, g 3 ) at z 3 (at time f 3 ). Denote by 
F y z the corresponding probability measure. It is easy to see that the map 

</> : C x D ^ [0,1], (y,z)^F ytZ {W 1 [0,T 1 ]n...nW p [0,T p ] / 0} 

is continuous and strictly positive, and independent of r by Brownian scaling. Hence 

c := inf 6(y, z) > 0. 

We infer that 

V[N x ,r\E Xtr nH t , r ] >c>0. 
Hence, combing our results, for sufficiently small r > 

F[E x , r n iV,, r |#,, r ] = F[E Xir j H Xtr ]¥[N x>r \ E x>r n > cr ^ 2 >-> 2 )+ £ , 

and this completes the proof as e > was arbitrary. o 
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3. Simulations 

To get hold of those exponents which we could not determine explicitly, we have performed Monte 
Carlo simulations. This has successfully generated conjectures in the p = 2 case, see Duplantier and 
Kwon [DK88j . Li and Sokal [LS90] and Burdzy, Lawler and Polaski [BLP89] . 

3.1 The general scheme. Before we list and analyse the simulated data, we explain how we got 
it. Fix positive integers p and n\, . . . ,n p . The aim is to get an estimate on ?(ni, . . . , n p ). Instead 
of Brownian motions we simulate two-dimensional symmetric nearest neighbour random walks. As 
it reduces computing effort, we work with boxes rather than with discs. (For comparison we have 
performed some of the simulations also with discs and there was no significant difference in the results.) 
First we fix an increasing sequence of box half-lengths Lq, . . . ,Lk (in most cases L k+ i = Ll-1 " 
and the maximal value m = Ll restricted to 20000, 40000 or 80000) and the sample size N of the 
simulation. 

Step 1. We start m + . . . + n p independent random walks at the origin € Z 2 and stop each of them 
when it hits d({— Lq, . . . , Lq} 2 ). This defines the starting positions of the random walks. 

Step 2. Assume we are at level k (after Step 1 we are at level k = 1). Independently run the random 
walks until they hit d({— L k , . . . , L k } 2 ). Separately, keep track of the set A^i C {— L k + 1, . . . , L k — l} 2 
of points that are visited by the ith package of random walks before hitting d({— Lf., . . . , L^} 2 ) (after 
Step 1). 

If Afc i fl . . . n Ak p = 0, then we say that we have survived level k and we enter level k + 1 (that is, 
we perform Step 2 again with k replaced by k + 1). Otherwise we stop this sample and start a new 
simulation in Step 1. 

By Nk we denote the number of samples that have survived level k. Clearly, iVo = N. We should have 

N k /N^(L k /L r^-^\ 

Hence in a double logarithmic plot of log(A^) against log(L&) the points should be on a line with slope 
—q(ni, . . . , n p ). Linear regression then gives an estimate for the exponent ?(ni, . . . , n p ). 

As it turns out that a line can be fitted well only for large values of L/%, we have neglected the small 
values of L k in order to get a reasonable estimate for ?(ni, . . . ,n p ). In Figure 13.4.11 below we plotted 
the data points used for the linear regression with solid circles, the other points with hollow circles. 

As can be seen from Figure 13.4.11 for £(1,1) this gives a pretty good estimate of the exact value |, 
even with a moderate computing effort of about 2000 hours CPU time. However, for £(1, 1, 1) the 
points tend to lie on a straight line only for large values of L k and thus require 

(i) a large maximal box size m = Lk and thus a big computer memory of size (2m + l) 2 bytes in 
order to keep track of the visited points, 

(ii) a large sample size iVo in order that Nk ~ No ■ (Lk/^o)"^ 1,1 ' 1 ' is big enough to obtain reliable 
data from the simulation. 

Since the CPU time we need for each sample grows with m, (i) and (ii) imply that we need huge 
amounts of CPU time. Furthermore, with huge sample sizes and box sizes, we run into the order of 
the cycle length of the common 48 bit linear congruence random number generators. 

The computations were performed on different computers, mainly on two parallel Linux clusters at 
the University of Mainz on Opteron 2218 processors with 2.6GHz and on Opteron 244 processors with 
1.8GHz. The programme code is written in C. As random number generator we used drand64(), a 64 
bit linear congruence generator following the rule 

r n+1 = (ar n + c) mod 2 64 
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with 

a = 6364136223846793005 and c = 1 

(see [Kn05l ppl06-108]). 

The linear regression method does not give a quantitative estimate on the statistical error. In order 
to get such an error estimate we did the following. Having in mind that the systematic error is large 
for small box sizes, we choose a minimal box number k m { n £ {1, . . . , K — 1} and neglect the data from 
all smaller boxes. Furthermore, we pretend that the asymptotics for is exact for k > k m i n , that is, 

(3.1) p Lk = C for all k > k min 

for some C > 0. In particular, the conditional probability to have no multiple intersections before 
leaving B^ k+1 given there is no multiple intersection before leaving Bi k is 

PL k+1 ( L k \ ? _ f 



Pk ■■= = t =: Q k 

PL k \L k+ i 1 

Here the likelihood function for the observation 

( N k min ,N kmin +i, ...,N K )=n:= (n kmin ,n kinin+l , . . . ,n K ) 

is 

K-l 



L n (,) = C(n) \[ p, n '+i(l-p,) n, " n,+1 



(3.2) '" fcmin 
v ' K-l 

= C{n) [] q^il-mT 1 -^ 1 

1—h ■ 

v "'mm 

for some C(n) > 0. The log-likelihood function is 

K-l 

(3.3) Cn(s) = log C(n) + ^ (jn+iS log(^) + (nj - n i+1 ) log (l - q\ 

1—h ■ 

L "-111111 

The maximum likelihood estimator (MLE) £ is defined by 

(3.4) C n U) = sup Cn(q). 

?>0 



We compute the derivatives 

K-l K-l 

(3-5) £' n (s)= n i+i l °g(qi) - ( n i- n i+i)-j 



log 



and 



l-l- I ! 

i ^niin 6 "'mill 



(3-6) C n {$) = - 2^ W - n l+ i) _ 2 . 

2-fc ■ V 1 <?;j 

Clearly, £"(?) < 0, hence ? i— > £ n (?) is strictly concave and thus <f is the unique solution of 

(3.7) £'„(<?) = 0. 

Hence, for given data, the MLE can easily be computed numerically (we used a Newton approximation 
scheme) . 
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Denote by <f no the MLE for sample size hq. By standard theory for MLEs, 0?)n o eN is consistent and 
asymptotically normally distributed. In fact, by Corollary 6.2.1 of |Le83j . 

(3.8) <f no n ^^° ? stochastically. 

Furthermore, by |Le831 Corollary 6.2.3], (<f)n is asymptotically efficient (that is, optimal) and by 
[Le83t Theorem 6.2.3] (with A/q,i the standard normal distribution) 



(3.9) yf no/(?) (<fn — n ^S° jVo,i in distribution. 

Here 



^r 1 / \ n™Y„.^2„? 



(3.10) I(,) = -E[£^(0|iV = l]= PLfemin Y, I II Pm)(l-Pl) (1 "f ( _ g ^/ 

2 ^min \ ^ ^min 



is the Fisher information for one sample. As we do not know the true value of ? and since we do not 
know pl. , we replace 1(0 by 

/„(,) = 

n 

By the law of large numbers In{s) n ^?° almost surely, uniformly in in compact sets. Hence 
by (|3.8p . we have /^(c) "^^l /(<j) stochastically. Hence we use 

(3.11) a 2 := -1/^(0 
as an estimator for the variance of <f and obtain 

(3.12) 1^2 Mo, i in distribution. 

a 

Concluding, an asymptotic 95% confidence interval for q is given by 

(3.13) [?-25, ? + 2a]. 

We have performed the simulations for the exponents <j(l, 1) and <j(2, 2) as benchmark problems, and 
then did the simulations on a larger scale for 

?(1,1,1), ?(1,1,2), 1,1,1), ?(1,1,1,2). 



3.2 Two-level scheme. The simulations turn out to be very time-consuming, especially for the 
exponents with a larger numerical value. In order to get a more efficient scheme in this situation 
consider the following simplification of the simulation scheme presented above: 

Assume there are only three box sizes, Lq (about 30), L\ (about (10 000) and L2 = 1L\. Then (13. 7|) 
can be solved explicitly and the maximum likelihood estimator for q is 

a _ log(n 2 /ni) 
' ~ log(2) ' 

In order to reduce the variance of <f we have to increase Ni, that is the sample size no- However, 
since it takes much CPU time to obtain a sample that contributes to Ni, we may wish to use this 
very sample as the starting point for a number m of trials running from box size L\ to L<i- Assume 
that x among these m trials have survived until L2 (that is, have reached the boundary of the L2-box 
without producing a multiple intersection), thenps = ^- is an estimator for the conditional probability 
of producing no multiple intersection until leaving the L2-box for the given realisation S of the paths 
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of all walks in the Li-box. Now we can prescribe the number n = n\ of "master samples" and for 
i = 1, . . . , n let Xi be the corresponding number of surviving trials and write pi := Xi/m. Hence for 



V-= — = %s] 

PL! 



we get the unbiased estimator 



n ^— ' 

1=1 



1 

7, 5^ 

The unbiased estimator for the variance of p is 

„ 1 n 

°"p = -7 TT 51 (Pi - ^) 2 - 

F n(n — 1) ^ 

From p and tr 2 we obtain the estimators for q and the variance a 2 of <f 
(3.14) ^=-^4 and ^- ^ 



log(2) ' "(log(2)p)2- 

We have employed this scheme for the exponents with numerical values larger than 2, and we explain 
now why it is more efficient in these cases. 

The expected time planar random walk needs to go from d({— L, . . . , L} 2 ) to d({— L — 1, . . . , L + l} 2 ) 
is of order L. The probability that a given sample ever reaches d({— L — 1, . . . ,L+ l} 2 ) is of order . 
Hence (if we stop the simulation as soon as the first multiple intersection is detected) the expected 
CPU time for each sample until box size L\ is of order 

Li 

E ^- 

L=L 

For q > 2 this sum is of order 1, for q < 2, it is of order L 2- '' '. Now the probability that a sample 
reaches box size L\ without producing a multiple intersection is of order . Hence the expected 
CPU time needed for simulating a "master sample" is of order L 2V? . On the other hand, each of the 
trials started from the master sample needs an expected CPU time of order L 2 . Hence for q > 2 we 
can run m = L\~ 2 trials without increasing the CPU significantly. 

In order to make a good choice for m, compute the variance of p 

Var[p] = n _1 Var[p s ] + E[p s (l -ps)\ < n^Varfps] + E[ps]. 

mn mn 

The quantities Var[ps] and E[ps] « 2"^ can be estimated from a test simulation as well as the expected 

CPU time T\ to produce a master sample and the expected time T2 used for each subsequent trial. 

Now it is an optimisation problem for the total CPU time n(T\ +TOT2) versus the variance Var[p]. For 

some of the simulations we have done test runs and solved the optimisation problem. Here m = 1000 

turned out to be a reasonable choice that we have then used in all simulations. 

We have performed the simulations according to this scheme with Lq = 30, L\ = 10 000, L2 = 20 000 
and m = 1000 for the exponents 

?(1,3,3), ?(2,2,2), ?(2,2,3), ?(2,3,3), ?(2,2,2,2). 
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3.3 Numerical results. We present our estimated values £ together with a statistical error of 2a. 
For the systematic error it is hard to make a good judgement. From the graphical representation of 
the results (see below) it seems that for <j(l, 1,2) the systematic error is of a smaller order than the 
statistical error. For 1, 2) and 1, 1) it is presumably of the same order. Finally, for <j(l, 1, 1, 1) 
and, even worse for <j(l, 1, 1, 1, 1) we seem to systematically underestimate the values. It would require 
a lot larger L max to get more accurate results. For that reason we have not taken too much effort to 
reduce the statistical error. However, we give the results of the simulations just to provide an idea of 
the possible values. 



exponent 




2a 


rigorous 


-^min 


^max 


no/10 6 


CPU 
time/h 


c(M) 


1.2502 


0.001 


5/4 


1069 


20 000 


500 


2 064 


?(2,2) 


2.9188 


0.0033 


f| = 2.9167 


163 


20 000 


40 000 


1879 


1, 1) 


1.027 


0.005 


[1/2, 5/4] 


18 575 


80 000 


60 


8 262 




1.2503 


0.0011 


[1,5/4] 


1069 


80 000 


200 


5 858 


^(1, 1,1,1) 


0.877 


0.006 


[1/4,5/4] 


39 813 


80 000 


20 


18 262 


1,1,2) 


1.02 


0.004 


[1/2, 5/4] 


27194 


40 000 


200 


35 212 


^(1,1,1,1,1) 


0.74 


0.02 


[1/8,5/4] 


27194 


40 000 


0.74 


1147 



Table 1. Numerical results obtained from the first simulation scheme. 



exponent 


<f 


2a 


rigorous 


n 


CPU 












time/h 


?(1,3,3) 


2.688 


0.01 


[2, (13 + v/73)/8] 


18100 


61860 


?(2,2,2) 


2.786 


0.01 


[2, 35/12] 


16 000 


47943 


<2,2,3) 


2.937 


0.01 


[2, 35/12] 


23 000 


116 888 


?(2,3,3) 


3.767 


0.057 


[2, 35/12] 


1000 


179 543 


?(2,2,2,2) 


2.664 


0.01 


[2, 35/12] 


16 000 


63496 



Table 2. Numerical results obtained from the second simulation scheme. 
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3.4 Detailed Data. 



3.4.1. Exponent <j(l, 1). The exact value 1) = 5/4 is known. This simulation is used as a bench- 
mark test for our simulation. 





n k 




Lk 




30 


500000000 




113 


109366745 


33 


455164209 




124 


97714439 


36 


414185142 




136 


87320799 


39 


379373384 




149 


78109962 


42 


349383901 




163 


69978568 


46 


315390855 




179 


62384176 


50 


286840826 




196 


55800459 


55 


257075021 




215 


49786852 


60 


232385705 




236 


44382636 


66 


207870728 




259 


39563995 


72 


187620511 




284 


35298660 


79 


168084821 




312 


31418279 


86 


151902122 




343 


27932867 


94 


136553134 




377 


24837149 


103 


122326905 




414 


22109889 



Lk 


rik 


455 


19660552 


500 


17483797 


550 


15525080 


605 


13788917 


665 


12253892 


731 


10890052 


804 


9669275 


884 


8589857 


972 


7631215 


1069 


6776772 


1175 


6020939 


1292 


5347118 


1421 


4747333 


1563 


4214131 


1719 


3741150 



Lk 


n k 




L k 


n k 


1890 


3323382 




7881 


557957 


2079 


2950258 




8669 


495180 


2286 


2620862 




9535 


439662 


2514 


2327160 




10488 


389839 


2765 


2066024 




11536 


345918 


3041 


1834523 




12689 


307046 


3345 


1628901 




13957 


272420 


3679 


1446024 




15352 


241798 


4046 


1283655 




16887 


214746 


4450 


1140213 




18575 


190486 


4895 


1012659 




20000 


173506 


5384 


898680 








5922 


797641 








6514 


708293 








7165 


628813 









Values used for the fit: L k = 1069 . . . 20 000. CPU time 2064h. 
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q (1 ,1 ) n = 500 000 000, L max = 20 000 
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5(1,1) = 1.2502 ±0.001 



Figure 2. Simulation for ^(1,1). The co-ordinate shows k m i n , the ordinate shows the 
corresponding <f with error bars. The vertical line indicates k m [ n = 40 which we chose 
for our estimate of <f. The horizontal line shows the true value. 



BROWNIAN INTERSECTION EXPONENTS 15 



3.4.2. Exponent ?(2,2). The exact value c(2, 2) 
serves as a benchmark for our simulations. 



35/12 = 2.91666 ... is known. Also this simulation 



T 


n k 




T 

^k 


n k 


30 


40000000000 




113 


459313243 


33 


27956276949 




124 


347384944 


36 


19934507109 




136 


263528000 


39 


14769670878 




149 


200799711 


42 


11270896745 




163 


153819037 


46 


8156016609 




179 


116600065 


50 


6108280379 




196 


89210485 


55 


4423365460 




215 


67932955 


60 


3315611015 




236 


51656232 


66 


2432628801 




259 


39313221 


72 


1842369408 




284 


30007400 


79 


1375726309 




312 


22780638 


86 


1056545535 




343 


17265563 


94 


803537797 




377 


13094893 


103 


607993657 




414 


9961095 





n k 


455 


7559087 


500 


5737717 


550 


4343548 


605 


3288311 


665 


2496057 


731 


1893876 


804 


1434709 


884 


1087314 


972 


823685 


1069 


624023 


1175 


473832 


1292 


359121 


1421 


271557 


1563 


205432 


1719 


155585 





n k 






n k 


1890 


117893 




7881 


1872 


2079 


89442 




8669 


1450 


2286 


67757 




9535 


1108 


2514 


51314 




10488 


853 


2765 


38803 




11536 


650 


3041 


29341 




12689 


479 


3345 


22363 




13957 


348 


3679 


16949 




15352 


266 


4046 


12813 




16887 


193 


4450 


9738 




18575 


151 


4895 


7339 




20000 


123 


5384 


5571 








5922 


4218 








6514 


3229 








7165 


2477 









Values used for the fit: L k = 605 . . . 20000. CPU time 1879h. 



S (2,2) n = 40 000 000 000, L max = 20 000 
J l l L 




<;(2,2) = 2.91 88 ±0.0032 



Figure 3. Simulation for c(2, 2). The co-ordinate shows fc m in, the ordinate shows the 
corresponding <f with error bars. The vertical line indicates k m [ n = 34 which we chose 
for our estimate of The horizontal line shows the true value. 
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3.4.3. Exponent <j(l, 1, 1). The exact value of <r(l, 1, 1) is unknown. 





n k 




Lk 


rik 


30 


60000000 




149 


25859377 


33 


59710616 




163 


24028791 


36 


58947709 




179 


22226328 


39 


57896946 




196 


20584233 


42 


56673833 




215 


19009323 


46 


54898618 




236 


17526913 


50 


53061195 




259 


16147690 


55 


50777396 




284 


14873454 


60 


48570208 




312 


13666336 


66 


46070175 




343 


12537025 


72 


43747356 




377 


11494795 


79 


41266693 




414 


10540910 


86 


39014779 




455 


9652748 


94 


36696389 




500 


8835893 


103 


34372986 




550 


8076259 


113 


32097711 




605 


7376857 


124 


29906035 




665 


6740503 


136 


27820941 




731 


6155608 



L k 


njt 


804 


5614962 


884 


5121770 


972 


4670981 


1069 


4256241 


1175 


3879722 


1292 


3534606 


1421 


3218775 


1563 


2930010 


1719 


2666485 


1890 


2426899 


2079 


2208165 


2286 


2009516 


2514 


1827541 


2765 


1661614 


3041 


1510468 


3345 


1372149 


3679 


1246493 


4046 


1132343 



L k 


ri k 




Lk 


rik 


4450 


1028587 




24722 


179046 


4895 


934379 




27194 


162421 


5384 


848491 




29913 


147273 


5922 


770449 




32904 


133514 


6514 


699523 




36194 


121333 


7165 


635064 




39813 


109856 


7881 


576119 




43794 


99544 


8669 


523319 




48173 


90226 


9535 


474777 




52990 


81910 


10488 


430885 




58289 


74069 


11536 


391134 




64117 


67148 


12689 


354964 




70528 


60809 


13957 


321882 




77580 


54981 


15352 


291736 




80000 


53301 


16887 


264553 








18575 


239803 








20432 


217707 








22475 


197364 









Values used for the fit: L h = 18575 . . . 80 000. CPU time 8262h. 



CM 
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S(1,1,1) = 1.027±0.0048 



70 



80 



Figure 4. Simulation for <j(l,l,l). The co-ordinate shows fc m i n , the ordinate shows 
the corresponding <f with error bars. The vertical line indicates fc m i n = 70 which we 
chose for our estimate of <f. The horizontal line shows the estimated value. 
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3.4.4. Exponent 1, 2). The exact value of <r(l, 1, 2) is unknown. 



^k 


n k 




^k 


n k 


ou 


zuuuuuuuu 




149 


O^OoO I Go 


33 


lyoiooiyy 




163 


4yzluo / 


36 


193436538 




179 


44119702 


39 


187242650 




196 


39647836 


42 


180308394 




215 


35523866 


46 


170682110 




236 


31776899 


50 


161176218 




259 


28415711 


55 


149901066 




284 


25417521 


60 


139521065 




312 


22668175 


66 


128301060 




343 


20191522 


72 


118361321 




377 


17981958 


79 


108205442 




414 


16028002 


86 


99385254 




455 


14267282 


94 


90671582 




500 


12694594 


103 


82304628 




550 


11279842 


113 


74445095 




605 


10020648 


124 


67187208 




665 


8909164 


136 


60566796 




731 


7917614 





n k 


QUI 


l uzyyoo 




0^40000 


972 


5545792 


1069 


4923405 


1175 


4374033 


1292 


3885012 


1421 


3449618 


1563 


3062025 


1719 


2718548 


1890 


2415286 


2079 


2144282 


2286 


1904636 


2514 


1690316 


2765 


1499756 


3041 


1331441 


3345 


1181425 


3679 


1048523 


4046 


930691 



^k 


n k 




^k 


n k 


aa 

44DU 






24722 


y / zoo 


4895 


I oooo I 




27194 


86246 


5384 


651874 




29913 


76472 


5922 


578486 




32904 


67873 


6514 


513593 




36194 


60371 


7165 


456239 




39813 


53674 


7881 


405240 




43794 


47628 


8669 


359523 




48173 


42353 


9535 


319391 




52990 


37627 


10488 


283792 




58289 


33387 


11536 


251871 




64117 


29608 


12689 


223767 




70528 


26339 


13957 


198585 




77580 


23407 


15352 


176146 




80000 


22541 


16887 


156219 








18575 


138744 








20432 


123285 








22475 


109424 









Values used for the fit: L k = 1069 . . . 10 000. CPU time 5858h. 
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Figure 5. Simulation for 1, 2). The co-ordinate shows k m i n , the ordinate shows the 
corresponding q with error bars. The vertical line indicates k m { n = 40 which we chose for 
our estimate of <f. The horizontal line shows the conjectured value ?(1, 1, 2) = ^(1, 1) = 
5/4. 
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3.4.5. Exponent cr(l, 1, 1, 2). The exact value of 1, 1, 2) is unknown. 
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179 


98203442 


46 


194683952 




196 
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215 
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236 
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60 
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259 
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66 
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284 
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72 
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312 
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79 


163642608 
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377 
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414 
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113 
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1175 
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1292 
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1421 
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Values used for the fit: L k = 27 194, . . . , 40000. CPU time 35 212h. 
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Figure 6. Simulation for <-(l, 1, 1, 2). The co-ordinate shows fc m i n , the ordinate shows 
the corresponding <f with error bars. The vertical line indicates fc m ; n = 74 which we 
chose for our estimate of The horizontal line shows the estimated value. 
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3.4.6. Exponent cr(l, 1, 1, 1). The exact value of 1, 1, 1) is unknown. 
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Values used for the fit: L k = 39 813 . . . 80 000. CPU time 18 262h. 
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Figure 7. Simulation for 1, 1, 1). The co-ordinate shows fe m i n , the ordinate shows 
the corresponding <f with error bars. The vertical line indicates fc m ; n = 78 which we 
chose for our estimate of <f. The horizontal line shows the estimated value. 
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3.4.7. Exponent q (1,1, 1,1,1) 
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Values used for the fit: L k = 27 194 . . . 40 000. CPU time 1147h. 
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Figure 8. Simulation for ^(1, 1,1,1,1). The co-ordinate shows /c m in, the ordinate 
shows the corresponding £ with error bars. The vertical line indicates k m i n = 74 which 
we chose for our estimate of The horizontal line shows the estimated value. 



3.4.8. Exponent <r(l, 3, 3). The exact value of ?(1, 3, 3) is unknown. As it turns out that <j(l, 3, 3) > 2, 
we have performed simulations according to our scheme 2. That is, we have generated n master 
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samples of random walk paths that reach the boundary of the Lx-box (here L\ = 10000). For each 
such master sample i we have run m = 1000 trials and have counted the fraction pi of trials where 
the paths reached the boundary of the L2-box (with Li = 2L\). As n = 18100 we cannot give the 
complete data set p\ , . . . , p n but rather give the empirical mean and the standard deviation of p 



p = 0.155202983425414, a v = 0.000536918044881792. 



From this we compute 



with standard deviation 



We conjecture that 



<?(1,3,3) = 2.6877718045551 



a = 0.00499094143436367. 



13 ~\- /73 

q(l, 3, 3) = 3) = — — — = 2.693000 . 



We conclude with a histogram of the values pf 



(1,3,3) Histogram of p, 18100 samples 
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Figure 9. Histogram of the values pi for 3, 3). 
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3.4.9. Exponent ?(2,2,2). The exact value of <j(2, 2,2) is unknown. We have performed a simulation 
with the second scheme with N = 16000, n = 1000, L x = 10000, L 2 = 20000. Mean and standard 
deviation are 

p = 0.1449495, 6-p = 0.000497221297799643. 

From this we compute 

£(2,2,2) = 2.78637773802317 

with standard deviation 

a = 0.00494888703003405. 

3.4.10. Exponent <j(2, 2, 3). The exact value of <j(2, 2,3) is unknown. We conjecture 

35 

?(2,2,3) =?(2,2) = — = 2.916666... 

We have performed a simulation with the second scheme with N = 23000, n = 1000, L\ = 10000, 
L2 = 20000. Mean and standard deviation are 

p = 0.130559, o-p = 0.000444444142417374. 

From this we compute 

<?(2, 2, 3) = 2.93722618256156 

with standard deviation 

a = 0.00491116935805033. 

3.4.11. Exponent ^(2,3,3). The exact value of ?(2,3, 3) is unknown. We conjecture 

?(2, 3, 3) = ?(2, 3) = 47 + 2 4^ = 3.738334113 . . . 

We have performed a simulation with the second scheme with TV = 1000, n = 1000, L\ = 10000, 
L2 = 20000. Mean and standard deviation are 

p = 0.073458, d v = 0.00144828442088002. 

From this we compute 

?(2,3,3) = 3.76693657262376 

with standard deviation 

a = 0.0284439101500224. 

This simulation was particularly time consuming (179 543h CPU time) as the actual value of ?(2, 3, 3) 
is rather large and it thus takes a tremendous amount of time to generate each master sample. 

3.4.12. Exponent ^(2,2,2,2). The exact value of ^(2,2,2,2) is unknown. We have performed a sim- 
ulation with the second scheme with TV = 16000, n = 1000, Li = 10000, L 2 = 20000. Mean and 
standard deviation are 

p = 0.157732125, a v = 0.000521232849038418. 

From this we compute 

£(2, 2, 2, 2) = 2.66445157389522 

with standard deviation 

a = 0.00476745017196814. 
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